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ABSTRACT 


In  this  paper  a  unified  presentation  is  made  on  the  results  of  var- 
ious investigators  on  the  properties  of  random  communication  networks. 
These  results  are  interpreted  in  such  a  way  that  the  properties  may  be 
determined  by  using  a  digital  computer  with  the  application  of  the 
Monte  Carlo  method.  The  computer  program  is  written  and  tested.  Re- 
sults for  some  networks  are  compared  with  theoretical  values. 
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I.   RANDOM  COMMUNICATION  NETWORK 

In  this  paper  the  results  of  various  investigators  on  the  properties 
of  random  communication  networks  are  reviewed.   From  these  results  a 
mathematical  model  is  represented  for  a  general  random  communication 
network . 

A.   THE  CONCEPT  OF  A  RANDOM  COMMUNICATION  NETWORK 

A  communication  network  is  an  aggregate  of  message  centers  that  at- 
tempt to  transfer  information  to  one  another  over  a  wide  variety  of 
channels.  The  message  centers  could  be  fixed  as  the  headquarters  of  a 
regional  military  command,  mobile  as  ships,  aircrafts,  or  satellites. 
The  channels  could  be  a  radio  link,  a  telephone  line,  or  highway.   In 
normal  conditions,  the  topology  of  the  network  is  known  with  certainty: 
the  positions  of  ships  at  sea,  the  links  between  the  headquarters  with 
different  units.  Each  center  knows  its  relative  position  with  respect 
to  the  other  centers  exactly.   If  a  message  is  sent  from  a  central 
command,  the  authorities  know  exactly  how  many  and  which  units  will  re- 
ceive it.  The  network,  in  this  case,  is  called  a  deterministic 
communication  network.   But  during  the  hostilities,  some  units  could  be 
destroyed  by  the  enemy  without  even  having  the  chance  to  send  the  last 
word,  some  radio  link  could  be  jammed  by  the  enemy,  some  telephone  lines 
or  some  roads  could  be  damaged  without  having  any  means  to  determine  the 
extent  of  the  destruction.   Back  at  the  headquarters,  the  authorities 
are  faced  with  the  uncertainty  of  the  existence  of  the  different  units 
and  channels,  the  effectiveness  of  an  order  sent  to  the  battlefield  is 
not  as  high  as  before.   It  is  only  probable  that  such  unit  gets  the 
message.  The  communication  network  becomes  a  highly  probabilistic, 

random  netwrok  in  the  event  of  a  war. 
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One  finds  the  same  kind  of  randomness  in  many  other  branches  of 
research;   for  example,  the  contagion  of  a  disease  among  a  population, 
an  hereditary  defect  in  a  species,  the  spread  of  a  stimulus  in  the 
nervous  system,  or  the  spread  of  information  in  society,  to  name  a  few. 
All  these  phenomena  are  very  similar  to  the  random  communication  net- 
work in  that  they  have : 

a.  an  originator  which  originates  a  specific  message 

b.  several  recipients  which  convey  the  received  'message'  to  the 
others  via  a  medium  (atmosphere,  voice,  physical  contact,  etc.) 

c.  the  uncertainty  of  contact  between  two  individuals. 

All  these  phenomena  could  be  then  studied  under  the  same  topic  of 
the  random  communication  network.   Researches  in  these  phenomena  give 
much  insight  in  the  problem  of  propagation  message  through  a  random 
network. 

As  it  is  said,  the  communication  between  two  message  centers  A  and 
B  could  be  assumed  only  with  a  certain  probability.  The  latter  measure 
is  a  function  of  many  parameters  which  are  for  example,  the  time  t 
when  the  message  is  sent,  the  time  T.   since  the  message  is  received  at 
station  A,  the  distance  d  between  stations,  the  electronic  equipment  of 
each  station,  to  cite  a  few.   Therefore  one  could  write,  in  general: 

Prob  (contact  A  to  B)  =  p  (A,  B,  t,"C  ,  d,  .  .  .). 

Due  to  the  complexity  of  the  problem,  without  reasonable  assumptions 
about  the  parameters  to  simplify  the  problem,  one  could  not  solve  it. 
Following  is  the  brief  survey  of  the  contributions  of  various  workers. 
Each  of  them  has  solved  partially  the  problem  with  different  assumptions 


about  the  parameters  which  determine  the  probability  of  contact  between 
two  stations,  two  individuals,  whatever  is  the  case. 

B.   THEORETICAL  DEVELOPMENT 

Bailey  [1]  studied  the  spread  of  an  epidemic  in  an  homogeneous 
population  with  size  n  where  the  probability  of  one  new  infection  tak- 
ing place  in  the  interval  dt  is  given  by 

Prob  (one  infection  in  dt)  =  y  (n-y+l)dt, 

where  y  =  the  number  of  susceptibles. 

From  the  above  equation,  the  mean  m(t)  of  the  number  of  infected  people 

and  the  epidemic  curve,  (i.e.,  the  rate  of  change  with  respect  to  time 

of  the  mean  of  the  number  of  infected  people)  dm(t)  are  obtained  to  be 

dt 


m(t)  =  e-10t(810t  -  234)  +  e-18  ^4410  t  -  902)  +  .  . 


z  =  dm£t)_    =  e-l0t(-8100t  _  3156-j  +  e-18t(79,380t  -  20650)  + 
dt 


Landau  and  Rapoport  [2]  assumed  the  probability  of  conta  ct  between 
any  two  individuals  is  the  same  for  any  pair  and  the  probability  of 
transmission  of  the  disease  depends  both  on  the  time  t  of  the  whole 
process  and  also  ont,  the  time  since  the  particular  affected  indivi- 
dual acquired  the  disease.  Thus, 

Prob  (one  infection)  =  p(t,X) 
The  differential  equations  giving  the  number  of  individuals  who  become 
affected  and  the  rate  are: 


dx(t)  =  <*(N  -  x) 
dt 


(t't)  +  Jo      dxU)  Vit''t   -A)d/\]  C1-1) 


X0P 


dzCt)  =  1  dx(t)  =  *N(1  -  z) 


^L°_=P(t,L)   (^  cU  p(t,t-A)  <uj 


dtNdt  LN  ^od 

(1.2) 

where  N  =  number  of  individuals  in  the  population 

3Q[t)  =  total  of  individuals  who  have  become  affected  up  to  time  t 
X0   =  number  of  individuals  who  become  affected  at  t  =  o,  the 

initial  time  of  the  process 
c(  =frequency  of  contacts,  which  is  assumed  constant.  .  More  pre- 
cisely ot,X(t)  [N  -  X(t)]  is  the  number  of  contacts  per  unit  time 
between  affected  and  unaffected  individuals. 

This  model  is  more  realistic  than  that  of  Bailey  because  it  is  true 
that  in  the  case  of  a  disease  epidemic,  the  infectiousness  of  the  path- 
ogenic virus  may  increase  or  decrease  with  the  time  of  the  process  and 
the  time  since  it  infects  the  body  of  the  individual. 

The  equations  (1.1)  and  (1.2)  are  difficult  to  solve  in  the  general 
case,  therefore  the  investigators  made  some  assumptions  and  solve  them 
for  the  following  cases: 

a.  The  probability  p(t,x)  is  a  function  of  t  alone 

b.  The  probability  p(t,r)  equals  e~^c 

c.  The  probability  p(t,x)  is  function  oft.  only 

d.  The  probability  p(t)  is  constant  in  a  finite  interval 

e.  There  is  a  lag  in  the  transmission  of  contagion. 

Landahl  [3]  treated  the  spread  of  some  information  as  a  flow  of 
'particles'  which  execute  random  motions  over  a  population  of  individuals 
and  which  may  multiply  or  disappear.   Equations  are  derived  for  the 
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density  of  these  particles,  and  the  number  of  individuals  through 

> 

which  the  'particles'  have  passed,  is  calculated.  The  results  are  ap- 
plied to  the  following  cases: 

a.  Uniform  spatial  distribution  with  multiplication  factor  decrea- 
sing with  time  because  of  loss  of  interest. 

b.  Multiplication  factor  is  constant,  but  the  rate  of  spread  de- 
creases with  multiple  hearings. 

c.  One  dimensional  region  with  a  small  starting  region  with  or 
without  an  absorbing  barrier.  An  absorbing  barrier  corresponds  to  the 
case  where  the  individuals  listen  to  the  information  but  do  not  repeat 
it. 

d.  Two-dimensional  region  with  absorbing  barrier. 

e.  Continuous  sources  of  information  within  a  small  region  in  one 
dimension. 

f.  Uniform  spatial  distribution  in  which  individuals  do  not 
respond  to  more  than  one  hearing. 

Prihar  [4]  studied  a  mobile  communication  network  where  the  trans- 
mitting and  receiving  stations  are  in  motion.  He  used  the  classical 
Kinetic  Theory  of  gases  to  compute  the  expected  total  number  of  vehi- 
cles contacted  per  unit  time  in  two  cases: 

a.  The  search  for  contact  is  continuous  and  the  antenna  is  omni- 
directional . 

b.  The  search  is  intermittent  and  the  antenna  is  directional. 
Finally,  Mattei  [5]  assumed  the  probability  of  contact  between  any 

pair  of  stations  is  constant  throughout  the  network. 

Prob  (contact  between  any  pair)  =  p  =  constant. 
Several  formulas  were  derived  in  closed  forms. 


The  mean  number  of  contacted  nodes  at  the  1st  step  is 

Z(l)  =  (n  -  l)p. 
The  mean  number  of  contacted  nodes  at  the  2nd  step  is 

Z(2)  =  (n  -  l)q  [1  -  (1  -  p2)n-2]. 
The  marginal  probability  of  the  number  of  new  nodes  at  each  step  is 


y(0)     yOJ        y(k-2) 
Prob[Z(k)  =  Z(k)]  =  5Z       J"       •••  > 


z(l)  =  0  z(2)  =  0      z(k-l)  =  0 

j(zuj)    [i-q«»)]«Ci)    [q»tP)]y(i)} 

(yzk(k)M  [i  -  ,«(k-«]«W  [qKk-i)]y00) 

k 
y(k)  =  n  -J~   z(i)      q  =  1-p  (1.3) 

where  n  =  number  of  nodes  in  the  network 
k  =  order  of  step 

C.   THE  PROPOSED  MATHEMATICAL  MODEL 

Although  with  the  model  of  Mattei,  one  could  derive  many  quantities 
of  the  network  in  closed  and  neat  forms,  the  model  is  not  realistic  in 
real  life  as  pointed  out  in  the  beginning  of  this  paper.  The  probabil- 
ity of  contact  between  any  pair  of  stations  is  far  from  constant.   It 
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depends  on  their  distances  to  the  center  of  attack,  the  size  of  the 
electronic  equipments,  the  geographical  positions.   In  the  event  of  a 
war,  the  capital,  the  missile  silos,  the  industrial  centers  are  more 
vulnerable  than  a  submarine  under  the  sea  or  a  communication  satellite 
in  orbit. 

In  this  paper,  a  model,  in  which  the  probability  of  contact  between 
any  pair  of  stations  could  be  varied  at  will,  could  be  assigned  to  each 
link,  once  all  circumstances  have  been  taken  into  account,  is  proposed 
and  solved.  The  solution  is  obtained  by  simulation  on  the  digital 
computer.  The  computer  program  allows  the  user  to  assign  the  probabil- 
ity of  contact  between  all  pairs  of  stations. 
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II.   SIMULATION  ON  A  DIGITAL  COMPUTER 

The  aim  of  this  chapter  is  to  obtain,  by  simulation  on 
the  digital  computer  using  the  Monte  Carlo  method,  the  de- 
sired quantities  of  a  general  random  communication  network, 
in  which  the  link  probabilities  could  be  all  distinct. 

A.   THE  MONTE  CARLO  METHOD 

The  Monte  Carlo  method  consists  of  solving  various 
problems  by  means  of  the  construction  of  some  random  process 
for  each  such  problem,  with  the  parameters  of  the  process 
equal  to  the  required  quantities  of  the  problem.   These 
quantities  are  then  determined  by  means  of  the  observations 
of  the  random  process  and  the  computation  of  its  statistical 
characteristics  which  are  approximately  equal  to  the  re- 
quired parameters. 

For  example,  the  required  quantity  x  might  be  the  math- 
ematical expectation  E(tt)  of  a  certain  random   .   The  Monte 
Carlo  method  for  determining  the  approximate  value  of  the 
quantity  x  consists  of  an  N-fold  sampling  of  the  value  of 
the  variable  V]  in  a  series  of  independent  tests,  (that  is, 
hi,  ^1  2 ,  ...   )  and  the  computation  of  their  mean  value 


1* 


n.  i  ^  t 


"»V 


N 
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If  the  number  of  tests  is  large,  one  has 

F^EC*]  )  =  x. 

It  is  appropriate  to  ask  the  question  why  one  does  not 
try  to  derive  an  analytical  expression  for  the  quantity  x 
from  the  properties  of  the  problem  instead  of  designing 
a  random  process  which  imitates  the  reality.   The  answer 
is  that  many  problems  in  the  physical  world  are  so  complex 
that  they  do  not  have  closed  formula  for  the  desired  quan- 
tities . 

Suppose,  for  example,   that  the  required  quantity  x 
is  the  definite  integral  of  a  function  f(y)  taken  over 
the  interval  (a,  b) ,  that  is 


■/ 


f(y)dy. 


m 


a- 


I; 


Fig.  2.1.   Illustrating  equation  (2.2) 
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If  the  expression  of  f(y)  is  not  simple  (or  even  worse, 
if  the  function  f(y)  does  not  have  a  closed  formula  but 
is  known  only  point  by  point) ,  x  could  not  be  evaluated 
analytically.   The  Monte  Carlo  method  may  be  used  as  will 
be  described  next. 

Suppose  we  have  a  device  called  a  random  number  gen- 
erator which  generates  two  random  variables  £  and  Y[ 
uniformly  distributed  ever  the  intervals  (a,  b)  and  (o,  c) , 
respectively.   The  Compound  probability  density  function 
is  given  by 


P^Vj   (y,z)  = 


c(b-a) 


0 


if  a<y<b  and  0<z<c 
otherwise 


The  pair  (b>)])  will  be  generated  N  times  and  each  time,  the 
condition 


f(§i)  .<  *]i 


(2.1) 


will  be  tested.   If  condition  (2.1)  holds,  the  point 
(£ i ,  hi)  in  Fig.  2.1  is  inside  or  on  the  boundary  of  the 
computed  area.   Let  N   be  the  number  of  times  when  this 
condition  holds  in  N  tests,  we  have 


■/: 


f(y)dy  = 


,.   .  lim  (N  i 
c(b-a)      \- ) 


(2.2) 
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In  order  to  obtain  good  results,  the  number  of  tests 
must  be  large,  say,  in  the  range  between  1,000  and  100,000. 
This  is  only  possible  with  the  aid  of  digital  computers. 

The  problem  of  a  random  communication  network  where 
the  probability  of  survival  of  each  link  is  assumed  constant 
has  been  solved  adequately  by  Mattei  [5].   If  one  changes 
this  assumption   and  lets  the  said  probability  vary  to  fit 
certain  conditions  in  reality,  the  problem  becomes  untrac- 
table,  then  the  only  feasible  solution  is  to  use  a 
numerical  method  such  as  the  Monte  Carlo  technique.   The 
required  quantities  are:  the  average  number  of  contracted 
nodes,  the  average  number  of  newly  contacted  nodes  at  each 
step,  the  terminal  reliability,  the  probability  distribution 
of  new  nodes  at  each  step.   The  method  for  obtaining  these 
quantities  will  be  described  in  the  following  section. 

B.   THE  FLOWCHART  OF  THE  SIMULATION  PROCESS 

Suppose  that  one  has  a  net^work  consisted  of  n  stations 
connected  by  a  certain  number  of  links,  it  is  a  determinis- 
tic network  in  the  sense  that,  if  a  message  is  sent  by  a 
station  to  adjacent  stations  which  relay  it  to  the  others, 
one  knows  exactly  how  many  and  what  stations  have  received 
it.   Suppose  a  disaster,  which  could  be  a  nuclear  attack, 
an  earthquake,  etc.,  happens  to  the  network.   Before  any 
investigation  of  the  damage  could  be  made,   one  does  not 
know  which  links  have  survived.   If  a  message  is  sent,  it 
is  difficult  to  acertain  how  many  and  which  stations  will 
get  the  message.   The  best  one  could  do  is  to  assign  to 
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each  link  a  certain  probability  of  survival.   Let  P.   be 
the  probability  of  communication  between  station  i  and  sta- 
tion j,  the  probabilities  between  all  pairs  of  stations  are 
best  described  by  an  n  x  n  matrix  as 


P 


(0 

P 


21 


P     P 
*11    12 


23 


P  ,   P  0 
nl    n2 


.  .  P 


.  .  P 


In 


2n 


.  .Pn,  n-1   0 


y1 


where   P   is   called   the  matrix   of   probabilities. 
Since   one   is   not   interested   in   the   case   where   a 
station   sends   back    to    iteself   a   received 
message,  all  diagonal  elements  of  P  are  zero.   From  a  set 
of  n  stations,  one  could  form  2  ^  "   .  networks.   Fig.  2.2 
shows  two  of  them  for  n  =  5. 


/ 


/ 


<s> 


/ 


& 


~Q 


Fig.  2.2.   Samples  of  random  network 
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If  the  number  of  stations  increases,  say  n  =  21  stations, 
the  population  of  possible  networks  is  2^  ~.  2  .  7  x 

19  6 

10    ,  which  is  large.   Suppose  that  one  investigates  the 
propagation   of  message  via  the  originator  which  is  the 
station  1.   Starting  with  this  station,  one  generates  a 
random  variable  r  uniformly  distributed  between  0  and  1, 
and  check  the  condition 

Q 
<§>-      -<3> 


step  0 
z(0)=l 


step  1 
z(l)=2 


step  2 
z(2)=4 


step  3 
z(3)=2 


z(4)=0 


Fig.  2.3.   A  tree  of  propagation 


r  x<  P 

*  » 


2' 


(2.3) 


If  this  condition  holds,  there  is  a  contact  between  station 
1  and  station  2  and  one  increases  the  random  variable  z(l), 
which  represents  the  number  of  newly  contacted  nodes  at 
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step  1,  by  one  unit.   One  generates  again  a  random  number  r, 
tests  condition  (2.3)  for  P     and  increases  z(l)  by  one 

1  |  o 

unit  if  condition  (2.3)  holds,  and  so  on  until  one  arrives 

at  the  last  station,  which  is  station  21.   One  repeats  the 
same  procedure  with  each  station  of  z(l)  relative  to  the 
remaining  network,  the  number  of  stations  contacted  at  this 
second  step  is  designated  z(2).   One  repeats  the  same  pro- 
cedure with  each  station  of  z(2).   The  procedure  stops  at 
the  step  when  one  could  not  contact  any  more  stations. 

As  an  example,  consider  the  network  in  Fig  2.3,  stations 
5  and  17  have  been  contacted  at  step  1  so  that  z(l)  =  2. 
The  stations  which  have  received  the  message  are  discarded 
from  the  network,  these  are  stations  1,  5,  17.   Next,  to 
each  station  contacted  at  step  1,  the  same  procedure  made 
for  station  1,  is  repeated.   In  the  example,  stations  2,  6, 
8  and  station  14  have  been  contacted  respectively  by  station 
5  and  station  17  at  step  2,  z(2)  =  4.   Stations  4  and  13 
have  been  contacted  by  station  6  at  step  3,  z(3)  =  2.   The 
propagation  of  the  message  is  over  at  step  4  because  stations 
4  and  13  cannot  carry  the  message  any  further. 

The  tree  in  Fig.  2.3  summarizes  the  process  in  various 
steps.   At  any  time  the  computer  needs  to  keep  the  infor- 
mation of  only  one  tree,  which  saves  a  lot  of  memory 
storage.   To  start  the  next  sampling,  one  just  has  to  erase 
all  information  of  the  first  sampling  run  except,  of  course, 
the  results  z(l),  z(2),  ....   The  general  procedure  for 
the  second  sampling  is  the  same  as  for  the  first. 
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After  K  samplings,  the  average  number  of  new  nodes  at 
each  step  is  given  by 


Z(l)  = 


Zx(l)  +  Z2(l)  +  Z3(l)  +  .  .  .  +  Zk(l) 
_ 


Z  (2)  +  Z  (2)  +  Z  (2)  +  .  .  .  +  Z  (2) 
Z(2)  =  -i 2 2 E 


Z,(KMAX)  +  Z  (KMAX)  +  .  .  .  +  Z     (KMAX) 
Z(KMAX)  =  ~ rr-=- 


where  KMAX  =  the  highest  order  of  steps  in  the  tree 
The  average  number  of  contacted  nodes  is 

X  =  1(1)  +  1(2)  +  .  .  .  +  Y(KMAX) . 


Zx(l)  +  Zj(2)  +  .  .  .  +  ZX(KMAX) 


+  Z2(l)  +  Z2(2)  +  .  .  .  +  Z2(KMAX) 


X  = 


+  Zk(1)  +  Zk  (2)  +  .  .  .  +  Zk(KMAX) 

_ 

X(l)  +  X(2)  +  .  .  .  +  X(KMAX) 
19 


where  X(i)  is  the  average  number  of  contacted  nodes  up  to 
step  i. 

The  Kth  terminal  reliability  is  commonly  defined  as 


r(K)  =  m. 

n-  1 


and  the  weak  connectivity  is  defined  as 


y    =  X(KMAX) 
I        n 


Definition .   In  this  paper  the  link  density  of  the  network 
is  defined  as  the  ratio  of  the  number  of  links  between  sta> 
tions  to  the  number  of  emitting  stations. 

During  the  process  of  simulation,  the  probability  dis- 
tribution of  new  nodes  at  each  step  is  also  calculated  and 
displayed  in  an  n  x  n  matrix  whose  rows  represent  the  num- 
ber of  steps  and  columns  represent  the  number  of  new  nodes 
as  f ol lows  : 

#  of  contacted  stations 
0    1    2    3    ...    n 


f 


1 


order  of  step 


J: 


Fig.  2.4.   Probability  distribution  of  the  number  of  con' 
tacted  stations  at  each  step. 
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A  flowchart  of  the  simulation  is  drawn  in  detail  in 
Fig.  2.5.   The  list  of  the  variables  used  in  the  computer 
program  is  described  as  follows. 


IRMAIN 


vector  whose  components  have  the  value  1  or  0. 

If  the  ith  component  is  1,  it  means  station  i 

has  yet  been  contacted,  0  means  the  station 
ith  has  been  contacted. 


NEW 


vector  whose  components  have  values  1  or  0.   If 

the  ith  component  is  1,  it  means  the  station  i 

has  just  been  contacted  during  the  preceeding 
step . 


ICOPY 


vector  with  1  or  0  as  components,  if  the  ith 
component  is  1,  it  means  the  station  i  has  been 
contacted  during  the  current  step.   Before  to 
go  to  the  next  step,  one  has  to  make;   NEW  = 
ICOPY. 


ZVEC 


vector  whose  ith  component  represents  the  num- 
ber of  newly  contacted  stations  at  step  (i  -  1) 


XVEC 


vector  whose  ith  component  represents  the  aver- 
age number  of  contacted  nodes  at  step  (i  -  1)  . 


GVEC 


vector  whose  ith  component  represents  the 
(i  -  l)th  terminal  reliability. 
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RELIAB   =    vector  whose  ith  component  represents  the. 
(i  -  1)  terminal  reliability. 

=    vector  whose  ith  component  represents  the  prob- 
ability distribution  of  new  nodes. 

RNGE     =    vector  whose  components  determine  the  x,  y 

ranges  of  the  plot  in  the  subroutine  UTPLOT  to 
plot  the  average  number  of  new  nodes. 

PRDIST   =    square  matrix  whose  (i,  j)  component  represents 
the  probability  that  the  number  of  contacted 
nodes  is  j  at  step  (i  -  1)  . 

PROB     =    square  matrix  whose  (i,  j)  component  represents 
the  probability  of  communication  from  station 
i  to  station  j  . 

RLIAB    =    square  matrix  whose  ith  row  represents  the  ter- 
minal reliability  for  the  ith  case  (different 
probabilities,  different  standard  deviations, 
etc . . . ) . 

XXVEC    =    rectangular  matrix  whose  ith  row  represents  the 
average  message  propagation  for  the  ith  case 
(different  probabilities,  different  standard 
deviations,  etc...). 
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KSTEP    =    integer  which  represents  the  order  of  the  step. 

KHECK    =    integer  whose  value  is  1  or  0  used  to  initialize 

KHECK1   =    ICOPY  to  (0)  or  to  compute  the  number  of  emitting 
stations . 

SUMND    =    number  of  contacted  nodes  up  to  step  K. 

ZNODES   =    number  of  contacted  nodes  during  step  K. 

KMAX     =    highest  order  of  step  attained  during  the  pro- 
cess of  simulation. 

KKMAX    =    KSTEP  +  2  used  in  the  computation  of  PRDIST 
minor  variable. 

CONECT   =    represents  the  weak  connectivity  of  the  network. 

NDEMIT   =    the  number  of  emitting  stations  during  one  sim- 
ulation . 

LABEL    =    used  in  the  subroutine  DRAW. 

MC       =    integer  used  in  the  subroutine  DRAW. 

IX       =    integer  used  in  the  subroutine  0VFL0W  and  RANDOM. 
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A       =    the  average  number  of  contacted  stations  at 

the  1st  step  in  the  model  of  Mattei  or  defined 
as  n/A,  standard  deviation  in  the  general  model. 

S        =    the  standard  deviation  of  destruction  in  the 
general  model,  defined  as:   S  =  n/A. 

X  and  Z  =    Normalized  variables,  defined  as  X  =  I/S  and 
Z  =  J/S. 

K        =    the  number  of  tests  in  the  process  of  simulation. 

P  and  PX  =   outputs  of  the  subroutine  NDTR  which  are  areas 
under  the  normal  curve. 

The  following  subroutines  of  the  NPGS  computer  center  library 
are  used  in  the  program. 

UTPLOT   =    to  plot  the  average  number  of  new  stations  curve. 

DRAW     =    to  plot  the  probability  distribution  of  con- 
nected stations,  the  terminal  reliability  and 
the  average  number  of  contacted  stations. 

OVFLOW,  RANDOM  =  used  to  get  a  random  variable  (0,1). 
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NDTR     =    to  obtain  the  area  under  the  normal  curve  as 
probability  of  communication  between  two  sta- 
tions . 

To  investigate  the  general  model  where  the  destruction 
is  assumed  normally  distributed  around  the  center  of  attack, 
block(2)  has  to  be  replaced  by  block(l). 
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SIMULATION  OF  A  RANDOM  COMMUNICATION  NETWORK 


IRMAIN  =  (11 
NEW  =  (00  . 
ICOPY  =  (00 
KSTEP  =  0} 
SUMND  =1.0, 
KKMAX  =  0, 
NDEMIT  =  0, 


.  .  1)   ZVEC  =  (10  ...  0) 
.  0)      XVEC  =  (10  ...  0) 
.  .  0)    GVEC  =  (0  ...  0) 
KHECK  =  0,     KHECK1^=  0 

ZNODES  =  0.0,     KMAX  =  0 
CONECT  =  0.0,     DNSITY  =  0.0 
MC  =1,     IX  =  102530469 
A=S=Z=X=P=PX=  0.0 

RNGE  =  (20.0,  0.0,  21.0,  0.0),  K  =  38,000 

1 


CALL  OVFLOW 


/" 


X 


N 


RLIAB  <r-  (0) 

GVEC  4~  (0  1  2  .  .  N) 

XXVEC  <t-  (0) 


"(D 


Fig.  2.5.   The  Flowchart  of  the  simulation. 
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PLOT:  TERMINAL 
RELIABILITY  AND 
PROPAGATION 
CURVES 


BLOCK  (1) 


PLOT  PROBABILITY 

DISTRIBUTION 

CURVE 


\ 


BLOCK  (2) 


A  <-  A  t-  1.0   PROB  (I, J)  = 

S  c-21/A   PROB  (J, I) 

X  v-  I/S   PROB  (I, J)  =  0 

CALL  NDTR    (I  =  J) 

Z  <-  J/S 

CALL  NDTR 

PROB  (I, J)  <-  (PX  -  0.5) (P-  0.5) 


A  <.-  A  +  1.0 


PROB  <-  (Ji_) 
N-l 

PRDIST  <.-  (0) 


© 


PRDIST  (1,2)  =  M 
IRMAIN  (1)  =  0 
NEW  (1)  =  1 
ICOPY  (1)  =  1 


-  G 


Q 


•  ' 

loo 

M<-1 

m<;k 

© 

M  c-  M  +   1 

Fig.  2.5.   continued, 
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© 


DNSITY 

«-0.0 

CONECT 

«-  0.0 

ZVEC  «- 

(0) 

RELIAB 

<-(0) 

DNSITY 

MESSAGE,  SPEED 
PROPAGATION 
PRDIST 
ZVEC  ^  " 


CALCULATION  OF: 
TERMINAL  RELIABILITY, 
PROPAGATION  SPEED, 
MESSAGE  PROPAGATION, 
PR.  DISTRIBUTION 

DNSITY   DNSITY/K 

ZVEC   ZVEC/K 


Fig.  2.5.   continued. 
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© 


(J) 


© 


1^1 


I  <~I  +  1 


I  4  N 


LA 

-* (    NEW(I)  =  0 


KHECK1  <r-  0 


J   1 


J   J  =  1 


J   N 


< 


IRMAIN(J)  =  0 


j 


F 


CALL  RANDOM 


"^ 


RANDOM   PROB(I,J) 


KHECK1  =  1 


3 


NDEMIT   NDEMIT  +  1 
KHECK1   1 


/ 


tL 


KHECK  =  1 


;F 

ICOPY 

(0) 

KHECK 

1 

© 


Fig.  2.5.   continued. 
29 


© 


© 


IRMAIN(J)    0 
ICOPY(J)    1 
SUMND   SUMND  +1.0 
ZNODES   ZNODES  +1.0 


f 


IRMAIN  =  (0) 


J 


f  T 


KSTEP  v  KSTEP  +  1 


^ 


t 


KMAX  >  KSTEP 


vF 


D 


KMAX   KSTEP 


ZVEC (KSTEP  + 

1)*- 

ZVEC (KSTEP 

+ 

1) 

+ 

ZNODES 

PRDIST (KSTEP 

+  1. 

ZNODES  +  1 

•0) 

•'-.  - 

- 

PRDIST(KSTEP 

+  1, 

ZNODES  +  1 

0) 

+ 

1 

.0 

ZNODES,.  0.0 

KKMAX,.-  KSTEP  +  2 


Y 


16  /—KKMAX 


16/— 16  +  1 


16  $  N 


Y 


PRDIST (16,1)  <- 
PRDIST(I6,1)  +  1.0 


-/ 


Fig.  2.5.   continued 


&' 
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a 


I© 


I 


NDEMIT  =  0 


\  v 


DNSITY 

(SUMND  -  1.0) 
NDEMIT 

NDEMIT   0 


=  DNSITY 


KSTEP   0 

KHECK   0 

SUMND   1.0 

IRMAIN    (1) 

NEW   (0) 

ICOPY   (0) 

4.-.. .   - 

<  "■  ■ 

Fig.  2.5.   continued. 
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© 


->— 


^Z<—\      NEW  t    I  COPY 


NEW  e-  I  COPY 

KSTEP  ^  -  KSTEP  -  1 

i 
: 

KMAX  <s  KSTEP 


KMAX   KSTEP 


KHECK  --  0 

ZVEC (KSTEP  +  1)  =  ZVEC (KSTEP  +  1)  +  ZNODES 

XVEC (KSTEP  +  1)  =  XVEC (KSTEP  +  1)  +  SUMND 

ZNODES^-  0 

. 

PRDIST(KSTEP  +  1,  ZNODES  +  1)  ^— PRDIST(KSTEP+1 , 

ZNODES  +  1) 

+  1.0 

Fig.  2.5.   continued 
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C.   COMPUTER  PROGRAM 

The  computer  program  is  written  to  analyze  a  network  of  21  sta- 
tions, however  this  number  could  be  increased  at  will.   First  one  has 
initialize  the  matrix  of  probabilities  of  communication  between  sta- 
tions to  desired  values.  The  probability  distribution  function  could 
be  discrete  or  continuous. 

For  example,  in  the  case  where  the  probability  of  communication 
between  any  pair  of  stations  is  constant,  it  suffices  to  write 

Table  2.1.   Fortran  statements  for  the  model  where  p  is  constant 


DO  16       I  =  1,  N 
DO  15      J  =  1,  N 
IF  (I.  EQ.  J)  GO  TO  14 
PROB  (I, J)  =  A/(N  -  1) 
GO  TO  15 

14  PROB  (I, J)  =  0.0 

15  CONTINUE 

16  CONTINUE 


If  one  wishes  to  analyze  a  network  where  the  destruction  is  normally 
distributed  around  the  center  of  a  nuclear  attack,  it  suffices  to 
replace  table  2.1  by  the  following  table 
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Table  2.2.   Fortran  statement  for  a  normal  model 


S  =  21/A 

DO  16   I  =  1,N 

X  =  I/S 

CALL  NDTR  (X,  PX,  D) 

DO  15   J  =  I,  N 

IF  (I.  EQ.  J)   GO  TO  14 

Z  =  J/S 

CALL  NDTR  (Z,  P,  D) 

• 

PROB  (I, J)  =  (PX  -  0.5)  . 

(P  - 

■  0.5) 

PROB  (J, I)  =  PROB  (I, J) 

14 

PROB  (I, J)  =  0.0 

15 

CONTINUE 

16 

CONTINUE 

The  computer  program  could  be  used  to  analyze  networks  of  any  size 
N  by  redimensioning  the  variable  names  which  are  vector  arrays,  matrix 
arrays  to  the  size  of  the  network  and  enter  the  data  card  according 
to  the  changes.  The  program  prints  out  the  link  density  and  the  weak 
connectivity  of  the  network,  as  well  as  the  curves  of  the  average 
number  of  contacted  stations,  the  probability  distribution  function  of 
the  number  of  contacted  stations  at  each  step  and  the  terminal  relia- 
bility of  the  network. 
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III.   COMPARATIVE  ANALYSIS 

For  the  purpose  of  comparison,  the  computer  program  has  been  used 
to  investigate  a  random  communication  network  having  21  stations 
whose  probability  of  contact  is  assumed  constant.  Three  cases  have 
been  considered  where  the  probability  is  equal  respectively  to  0.05, 
0.10,  and  0.15.  The  sampling  size  is  calculated  to  have  an  error  of 
0.01  and  the  results  are  compared  with  the  theoretical  values. 

A.   THE  CALCUATION  OF  THE  SAMPLING  SIZE 

Suppose  that  one  station  tries  to  send  the  message  to  all  other 
stations.  The  expected  number  of  contacted  stations  is  given  by 

P"  =  (n  -  l)p.  (3.1) 

The  variance  is  given  by 

CT  =  (n  -  l)pq 

where  n  =  the  size  of  the  network 

p  =  the  probability  of  contact 
q  =  1  -  p. 

Equation  (3,1)  becomes 


P- 


=  1 


in  the  case  where  n  =  21  and  p  =  0.05, 
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Let  Z(l)  designate  the  random  number  of  contacted  stations  at  the 
first  step  in  one  sample.   Since  the  population  is  infinite,  the  mean 
and  the  standard  deviation  of  the  sampling  distribution  of  means  are 
given  by 


Kar  z(d  =^ 


and   ^ZM)  =-£= 


where  N  =  the  sampling  size. 

For  large  values  of  N  (N  ^  30)  the  sampling  distribution  of  means 
is  approximately  a  normal  distribution  with  mean  U  and  standard  de- 
viation G~ Kf7      •  The  accuracy  of  the  approximation  improves  as  N 
gets  larger  (Central  limit  theorem) . 

Suppose  that  one  would  like  the  value  of  Z(l)  differ  from  the 
population  mean  Lt  by  0.01  with  a  confidence  level  of  95.45%,  one  has 
to  satisfy  the  inequality 

|  Z(l)-^|  ^  2(TZ(1)  ^  o.Ol  (3.2) 

2^1)pq     «   10-2 

or      N  ^  4.104   .   (n-l)pq  =  4.104  .  (21-1) (0.05) (0.95) 

then    N  ^  3.8  x  104 

The  program  has  been  run  through  3.8  x  104  iterations.  The  results 
of  the  simulation  are  compared  with  the  theoretical  values  in  the  fol- 
lowing section. 
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B.   COMPARISON  BETWEEN  SAMPLING  AND  THEORETICAL  VALUES 

a.  Average  Numbers  of  Contacted  Stations  at  the  First  and  Second 
Steps 

The  average  numbers  of  contacted  stations  are  given  by  [5] 


Z(l)  =  (n-l)p 


Z(2)  =  (n  l)q  jl  -  (l-p2)n"2}. 


They  are  calculated  for  three  values  of  the  probability  p  and  compared 
with  the  results  of  the  sampling  in  tables  3.1  and  3.2.  The  difference 
between  sampling  and  theoretical  results  is  small. 

Table  3.1.   Values  of  Z(l) . 


._::._    ._. .-,. 

Probability 

Theoretical 

Sampling 

Difference 

0.05 
0.10 

1 

0.9978 

0.0022 

2 

2.0007 

-0.0007 

0.15 

3 

2.9870 

0.0130 

Table  3.2.  Values  of  Z(2) 


Probability 

Theoreti 

I 

cal 

Sampling 

Difference 

0.05 

0.8825 

0.8708 
3.1334 

0.0117 

0.10 

3.1289 

0.0045 
0.0306 

0.15 

5.9676 

5.9370 
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The  sampling  size  has  been  calculated  for  the  value  of  probability 
equal  to  0.05.  Table  3.1  shows  that  the  difference  between  the  sampling 
and  theoretical  values  z"(l)  -  1   is  0.0022,  less  than  0.01  as  it  is 
expected. 

b.  Average  Total  Numbers  of  Contacted  Stations 

Fig.  3.1  gives  the  numbers  of  stations  which  will  be  eventually 
contacted  for  different  values  of  probability.  They  are  tabulated  in 
Table  3.3  along  with  the  corresponding  theoretical  values. 


Table  3.3.   Average  total  number  of  contacted  stations 


Probability 

i 
Theoretical       Sampling 

! 
Difference 

0.05 

4.50            4.40 

0.10 

- 

0.10 

14.00           14.07 

-0.07 
0.10 

0.15 

19.00           19.10 
1 

The  theoretical  [5]  and  sampling  curves  giving  the  average  numbers  of 
contacted  stations  up  to  step  K  are  similar.   Both  show  that  the 
propagation  of  message  is  over  between  steps  4  and  6. 
c.  Terminal  Reliability 

The  terminal  reliability  curves  in  [5]  and  in  Fig.  3.2  are 
very  similar.  They  show  that  the  highest  reliability  occurs  between 
steps  2  and  3  and  is  negligible  from  step  6.  The  terminal  reliabil- 
ity is  the  measure  of  the  speed  of  propagation  of  the  message. 
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Table  3.4.   Values  of  highest  terminal  reliability 


Probabil. 

ity  ! 

Theoretical 

Sampling 

Difference 

0.10 

| 

0.175 

0.175 
0.32 

0.0 

0.15 

0.34 

0.02 

d.   Probability  Distribution  of  the  Number  of  New  Stations 
Contacted  at  Each  Step 

Fig.  3.3,  3.4,  and  3.5  give  the  probability  distribution  of 
new  stations  contacted  at  each  step,  they  correspond  to  the  formula 
(1.3),  they  are  very  similar  to  the  theoretical  curves  in  [5], 

In  summary,  the  computer  program  has  yielded  results  almost 
identical  to  the  theoretical  values  if  the  sampling  size  is  large, 
between  4  x  10  and  10  .   Since  the  program  was  written  without  any 
constraint,  any  assumption  on  the  probability  of  communication  between 
pairs  of  stations,  it  could  be  used  with  confidence  to  study  any  ran- 
dom network  on  the  condition  that  the  matrix  of  probabilities  is 
initialized  to  values  which  reflect  each  particular  circumstance.  As 
an  example,  in  the  next  chapter,  the  program  is  used  to  investigate 
a  communication  network  whose  destruction  is  supposed  normally  dis- 
tributed around  a  nuclear  explosion  or  an  earthquake. 
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Fig.  3.1.  Average  message  propagation 

x  scale  =  order  of  step  (2  units  per  inch) 

y  scale  =  number  of  contacted  stations  (3  units  per  inch) 
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Fig.  3.2.  Terminal  Reliability 

x  scale  =  order  of  step  (2  units  per  inch) 

y  scale  =  terminal  reliability  (0.1  unit  per  inch) 
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Fig.  3.3.   Probability  distribution  of  the  number  of  contacted 
stations  at  each  step. 

Probability  =0.05 

x  scale  =  number  of  contacted  stations  (3  units  per  inch) 

y  scale  =  probability  distribution  (0.1  unit  per  inch) 

k      =  order  of  step 
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Fig.  3.4.   Probability  distribution  of  the  number  of  contacted 
stations  at  each  step. 

Probability  =0.1 

x  scale  =  number  of  contacted  stations  (3  units  per  inch) 

y  scale  =  probability  distribution  (0.1  unit  per  inch) 

k      =  order  of  step 


49 


\  y^^C2 


a:3 


:c5 


Oi2 


Fig.  3.5.   Probability  distribution  of  the  numbers  of  contacted 
stations  at  each  step. 

Probability  =0.15 

x  scale  =  number  of  contacted  stations  (3  units  per  inch) 

y  scale  =  probability  distribution  (0.1  unit  per  inch) 

k      =  order  of  step 
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IV.   INVESTIGATION  OF  A  GENERAL  MODEL 

A.   REPRESENTATION  OF  THE  MODEL 

Consider  a  communication  network  of  n  stations  which  may  be  rep- 
resented by  an  n  x  n  matrix  whose  elements  are  1  or  0  as 


P. 


n 


12  3.. 
0  10  0  1 
0  0  10  0 
10  0  0. 
0 
1 


n 
1 

0 

1 


(4.1) 


Suppose  there  is  a  perturbation  near  station  1,  whose  effect  is  to 
destroy  the  links  of  communication  according  to  a  normal  distribution 
N(0,<r)  around  a  center.  This  means  that  the  destruction  is  maximum 
at  the  point  of  impact  and  decreases  radially  according  to  the 
equation 


P(x)  =f  _i_  e-V^2/<r2 
/o  ctVTtF 


dt. 
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Fig.  4.1.   Illustrating  the  equation  (4.2) 

where  p(x)  =  probability  of  survival  of  a  link  of  a  station 
x  =  distance  of  a  station  to  the  impact 
<7~  =  standard  deviation  of  the  destruction 

The  probability  of  communi cation  between  stations  i  and  j  is 
Pij  =  P(*i)  p(Xj) 
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2fl2 


.; i  _i_  \^  dt/^Ts e^  dt-      (4-2) 


0-V2^ 


Suppose  that  the  stations  are  numbered  according  to  increasing  dis- 
tance and  the  network  is  complete,  then  all  elements  of  the  matrix  P 
have  to  be  replaced  by  the  values  of  equation  (4.2).  The  value  of 
the  standard  deviation  measures  the  extent  of  the  destruction.  The 
simulation  was  made  for  three  cases  where  station  21  is  one,  two  and 
three  standard  deviations  away  from  the  impact. 

first  case  <j-  =  21  extensive  destruction 

21 
second  case  CT   =  —  =  10.5      medium  destruction 

third  case  /j-  =  ±L  =   7        minor  destruction 

3 

The  subscripts  of  the  stations  which  represent  their  distance  to  the 
impact,  have  been  expressed  in  standard  units  and  the  subroutine  NDTR 
has  been  used  to  initialize  the  matrix  P.  The  sampling  size  has  been 
taken  equal  to  3.8  x  10  .  The  results  are  summarized  in  the  following 
section. 

B.   RESULTS  OF  THE  SIMULATION 

When  the  farthest  station  is  one  standard  deviation  from  the  im- 
pact, i.e,  the  case  when  the  destruction  is  extensive,  the  message  is 
not  carried  to  many  stations.   On  the  average,  1.26  stations  receive 
the  message,  the  link  probabilities  being  very  low.   In  the  case  of 
medium  destruction,  the  farthest  station  two  standard  deviations  from 
the  impact,  4.17  stations  receive  the  message.  The  propagation 
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ceases  at  step  5.   7.45  stations  receive  the  message  when  station  21 
is  three  standard  deviations  away  from  the  impact.  The  above  results 
are  obtained  with  the  assumption  the  communication  network  forms  a 
complete  graph,  i.e.,  there  is  a  link  between  any  pair  of  stations. 
It  is  usually  not  the  case  in  real  life,  therefore  the  above  results 
are  only  the  upper  bound  of  number  of  stations  which  will  receive  the 
message  in  case  of  a  war. 

If  one  compares  with  the  model  where  the  link  probability  is  con- 
stant, one  finds  the  latter  model  always  gives  much  higher  results 
than  the  current  model  which  is  nearer  to  reality.  Table  4.1  summar- 
izes the  results  of  two  models. 


Table  4.1.  Number  of  contacted  stations  in  two  models 
compared. 


Probability 


Mattei 


0.05 
0.10 
0.15 


4.40 


14.07 


19.10 


Standard  Deviation    Normal  Model 


21 

10.5 
7 


1.25 

4.17 
7.45 


Fig.  4.2  shows  the  average  message  propagation  in  the  three  cases. 
The  weak  connectivity  of  the  normal  model  is  also  smaller  than  the 
Mattei  model,  it  is  equal  to  0.06,  0.20  and  0.35  when  the  farthest 
station  is  away  from  the  impact  one,  two  and  three  standard  deviations 
respectively.  The  link  densities  are  respectively  0.08,  0.36  and 
0.67.  All  results  point  out  that  in  the  normal  model,  the  propagation 
stops  rapidly.  The  terminal  reliability  in  Fig.  4.3  is  a  monotoni- 
cally  decreasing  function  in  the  first  case,  attains  the  maximum 
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value  at  step  3  for  the  second  and  third  cases.   It  means  that  the 
propagation  increases,  attains  the  culmination  then  decreases  and 
stops  in  average  at  step  5.   The  curves  of  probability  distribution 
function  in  Figs.  4.4,  4.5  and  4.6  of  the  number  of  contacted  stations 
at  each  step  show  also  small  values  compared  with  the  Mattei  model. 


55 


LI 
.-1 

t  i 
-1 

II 
£  J 

■ 

^     1 

IP 

( ; 
rj 

«r^  10.5 

c> 
ci 

<r-  24 

LI 

r> 
o 

JUw 


Q2<! 


a:s 


Fig.  4.2.  Average  message  propagation 

x  scale  =  order  of  step  (2  units  per  inch) 

y  scale  =  number  of  contacted  stations  (3  units  per  inch) 
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Fig.  4.3.  Terminal  reliability 

x  scale  =  order  of  step  (1  unit  per  inch) 

y  scale  =  terminal  reliability  (0.1  unit  per  inch) 
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Fig.  4.4.  Probability  distribution  of  the  number  of  contacted 
stations  at  each  step 

Standard  deviation  =  21 

x  scale  =  number  of  contacted  stations  (3  units  per  inch) 
y  scale  =  probability  distribution  (0.2  unit  per  inch) 
k      =  order  of  step 
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Fig.  4.5.  Probability  distribution  of  the  number  of  contacted 
stations  at  each  step 

Standard  deviation  =  10.5 

x  scale  =  number  of  contacted  stations  (3  units  per  inch) 
y  scale  =  probability  distribution  (0.2  unit  per  inch) 
k      =  order  of  step 
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Fig.  4.6.  Probability  distribution  of  the  number  of  contacted 
stations  at  each  step 

Standard  deviation  =  7 

x  scale  =  number  of  contacted  stations  (3  units  per  inch) 

y  scale  =  probability  distribution  (0.2  unit  per  inch) 

k      -  order  of  step 
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V.   CONCLUSION 

The  analysis  of  the  message  propagation  through  a  random  communi- 
cation network  has  been  studied  by  various  investigators.  To  be  able 
to  solve  the  problem,  each  has  adopted  different  assumptions  on  the 
probability  of  communication  between  stations.  The  main  effort  of 
this  paper  has  been  to  solve  the  same  problem  without  having  to  make 
any  simplifying  assumption  on  the  link  probability.  The  latter  could 
vary  to  fit  a  particular  circumstance. 

A  flowchart  has  been  drawn  and  a  program  for  its  implementation 
on  the  digital  computer  has  been  written.  The  program  gives  the  de- 
sired quantities  of  a  random  communication  network  whose  link 
probabilities  could  be  varied  to  fit  a  particular  situation  by  the 
Monte  Carlo  technique.   The  program  has  been  tested  with  the  model 
whose  link  probabilities  were  assumed  constant  and  whose  characteris- 
tics could  be  derived  analytically.  The  results  were  found  in 
agreement  with  the  analytical  results.  The  program  has  also  been 
used  on  a  model  where  the  destruction  was  supposed  to  be  normally 
distributed  around  the  point  of  attack.  The  results  also  agree  with 
the  intuitive  idea  that,  as  the  effectiveness  of  the  destruction  of 
the  network  increases,  the  propagation  of  message  stops  rapidly. 

In  the  computer  program,  the  time  taken  by  the  message  to  travel 
between  any  pair  of  stations  is  assumed  constant,  which  is  not  true 
in  real  life.   Some  stations  take  longer  time  to  process  and  send  the 
information  than  the  others.   It  is  proposed  for  future  work  that  the 
time  t^ j ,  time  taken  by  station  i  to  process  information  and  send  it 
to  station  j,  could  be  varied  at  will.  With  this  second 
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generalization  added,  the  simulation  will  give  the  average  time  of  the 
propagation  of  the  message  and  what  routing  procedure  will  give  the 
fastest  communication  between  various  stations. 
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